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The boson images of fermion SO(2N+l) Lie operators have been given together with 
those of SO(2N + 2) ones. The SO(2N+l) Lie operators are generators of rotation 
in the (2iV+l)-dimensional Euclidian space (N: number of single-particle states of 
the fermions). The images of fermion annihilation-creation operators must satisfy the 
canonical anti-commutation relations, when they operate on a spinor subspace. In the 
regular representation space we use a boson Hamiltonian with Lagrange multipliers to 
select out the spinor subspace. Based on these facts, a new description of a fermionic 
50(2iV + 2) top is proposed. From the Heisenberg equations of motions for the bo- 
son operators, we get the SO(2N + 1) self-consistent field (SCF) Hartree-Bogoliubov 
(HB) equation for the classical stationary motion of the fermion top. Decomposing 
an SO(2N+l) matrix into matrices describing paired and unpaired modes of fermions, 
we obtain a new form of the SO(2N+l) SCF equation with respect to the paired- 
mode amplitudes. To demonstrate the effectiveness of the new description based on 
the bosonization theory, the extended HB eigenvalue equation is applied to a super- 
conducting toy-model which consists of a particle-hole plus BCS type interaction. It is 
solved to reach an interesting and exciting solution which is not found in the traditional 
HB eigenvalue equation, due to the unpaired-mode effects. To complete the new de- 
scription, the Lagrange multipliers must be determined in the classical limit. For this 
aim a quasi anti-commutation-relation approximation is proposed. Only if a certain 
relation between an SO(2N-\i) parameter z and the N is satisfied, unknown parameters 
k and I in the Lagrange multipliers can be determined withuout any inconcistency. 
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1 Introduction 



The time dependent Hartree-Bogoliubov (TDHB) theory is the leading standard approx- 
imation in the many-body theoretical description of a superconducting fermion system pQ- 
[3]. The HB wave function (WF) for the fermion system represents Bose condensate states 
of fermion pairs. It is a good approximation for the ground state of the fermion system 
with a short-range pairing interaction that produces two-body correlations which is taken 
into account as a spontaneous Bose condensation of the fermion pairs. The fermion number- 
nonconservation of the HB WF is a consequence of the spontaneous Bose condensation 
of fermion pairs which causes a coherence in phases of superconducting (Bose condensed) 
fermions. The SO(2N) Lie algebra of the fermion pair operators contains the U(N) Lie 
algebra as a subalgebra. Here the SO(2N) and the U(N) denote the special orthogonal 
group of 2N dimensions and the unitary group of N dimensions, respectively (N: number 
of single-particle states of the fermions). The canonical transformation of the fermion op- 
erators generated by the Lie operators in the SO(2N) Lie algebra induces the generalized 
Bogoliubov transformation for the fermions [I] to vanish the dangerous term for arbitrary 
pairs. The TDHB equation is derived from the classical Euler-Lagrange equation of motion 
for the coset variables [5]. Usually the solutions of this equation provides the ground 

state of an even fermion system. For the odd fermion system, we must consider a one quasi- 
particle state on the HB WF in which paired and unpaired states are not treated in an equal 
manner. In this sense for such a system we have no TD self-consistent field (SCF) theory 
with the same power for the mean field approximation (MFA) as the TDHB theory. 

One of the most challenging problems in current studies of condensed matter physics is 
to obtain a theory suitable for the description of collective motions with large amplitudes 
in fermion systems with strong collective correlations. For providing a general microscopic 
framework for a unified SCF description for Bose- and Fermi- type collective excitations in 
those systems, a many-body theory has been proposed by Fukutome, Yamamura and one of 
the present authors (S.N.) based on the SO(2N+l) Lie algebra of the fermion operators [6]. 
An induced representation of an SO(2N+l) group has been obtained from a group extension 
of the SO(2N) Bogoliubov transformation for fermions to the SO (2 N+l) transformation 
group. We start with the fact that the set of the fermion operators consisting of the creation- 
annihilation and the pair operators forms a larger Lie algebra, the SO(2N+l) Lie algebra. 
The fermion Lie operators, when operating on the integral representation of the SO(2N+l) 
WF, are mapped into the regular representation of the SO(2N+l) group and are represented 
by Bose operators. The Bose images of the fermion Lie operators are expressed by closed 
first order differential forms. The creation-annihilation operators themselves as well as the 
pair operators are given by the finite Schwinger type boson representation [7J [8]. 

Embedding the SO(2N+l) group into an SO(2N+2) group and using the boson images of 
,50(2^+2) Lie operators, we have developed an extended TDHB (ETDHB) theory [5] and 
extended supersymmetric a- model [TU| ITT] . Particularly in the ETDHB theory, which is ob- 
tained from the Heisenberg equation of motion for the boson operators, paired and unpaired 
modes are treated on an equal footing. A static extended Hartree-Bogoliubov (EHB) theory 
is derived easily from the ETDHB theory. The EHB theory applicable to both even and odd 
fermion systems is a SCF theory with the same power for the MFA as the usual HB theory for 
even fermion systems. Based on these facts, a new description of a fermionic 5*0(2A^+2) top 
is proposed. We start from the Hamiltonian of the fermion system which includes, however, 
the Lagrange multipliers needed to select the physical spinor subspace. The EHB equation 
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is written in terms of appropriate variables representing the paired and the unpaired modes. 
The EHB eigenvalue equation is solved by a method parallel to the two-step diagonalization 
method for the usual HB eigenvalue equation [12]. We obtain a new eigenvalue involving 
unpaired-mode effects in contrast to the usual HB theory, which is unable to describe the 
unpaired modes. Through the coordinate transformations for space fixed and body fixed 
coordinate frames, the fermion SO(2N+l) Lie operators are expressed in terms of the quasi- 
particle expectation values (c-number) of them and the quasiparticle SO(2N-\l) Lie operators 
(quantum mechanical fluctuations). To treat the quantum mechanical fluctuations, we make 
successive coordinate transformations. The fluctuating Hamiltonian with the Lagrange mul- 
tipliers is given up to the first order in the quasi-particle SO(2N+l) Lie operators. We found 
that Sawada's eigenmode method [131 [H] is not appropriate to obtain the excitation energy 
due to the Hamiltonian because it contains the unpaired-mode amplitudes in a particular 
form. However, with the help of the traditional Bogoliubov's approach [I], we were able to 
derive a very simple expression for the excitation energy due to the fluctuating Hamiltonian. 
The unknown parameters in the Lagrange multipliers terms must be determined. The anti- 
commutators of the fermion Lie operators given in the first order differentials satisfy exactly 
the anticommutation relations, when they operate on the SO(2N+l) HB WF. However, this 
fact plays no role to determine the parameters. A determination of the parameters is possible 
if we demand instead that expectation values of the anticommutators by an approximate 
SO(2N+l) HB WF satisfy the anticommutation relations, i.e., the quasi anticommutation 
relation approximation for the fermions [9j[T5]. Under the approximation, the determination 
could be made successfully. This plays crucial roles for a unified self-consistent description 
of Bose-Fermi type collective excitations at regions very near z—\ (the case of non existence 
of unpaired modes) and z=0 (the case of largest contribution from unpaired modes) since 
they are the necessary and sufficient conditions to determine the unknown parameters in 
the Langrange multipliers. This means a great step to towards verification of the validity 
of the present new description. The above quantum fluctuations can also be studied in the 
framework of quantum group symmetry formalism PS]- [19]. Finally, we point out that it is 
an important and interesting problem to attempt a group theoretical approach to formation 
of the Lax pair of the SO(2N+2) top [H]-[22]. 

This paper is organized as follows. In Sec. 2, we recapitulate briefly the induced represen- 
tation of the SO(2N+l) canonical transformation group, the embedding of the SO(2N+l) 
group into an SO(2N+2) one and the introduction of the S yf^t? coset variables. In Sec. 3, 



we give a brief sketch of the derivation of the extended TDHB equation from the Heisenberg 
equation of motion for the boson operators. The extended TDHB theory is just the TD- 
SCF theory. In Sec. 4, a static extended HB equation is derived from the extended TDHB 
equation. The EHB eigenvalue equation is solved with the two-step diagonalization method. 
In Sec. 5, energies of classical motion and quantum mechanical fluctuation of the fermion 
SO(2NH) top are given. Sections 6 and 7 are devoted to a determination of the unknown pa- 
rameters in the Lagrange multipliers by the quasi- anticommutation relation approximation 
for the fermions. Finally, in the last section, we give some concluding remarks and further 
perspectives. In particular we will attempt a group theoretical approach to the formation of 
the Lax pair of the SO(2N-\2) top. In App. A, we give expressions for differential SO{2N J \2) 
Lie operators on the coset spaces. Further we provide useful formulas to evaluate the ex- 
pectation values of differential annihilation-creation operators. We characterize the action 
of the corresponding operators on the SO(2N+l) wave function. Throughout this paper, we 
use the summation convention over repeated indices unless otherwise stated. 
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2 Brief review of Bogoliubov transformation generated 
by SO(2N+l) Lie algebra of fermion operators and 
its embedding into SO(2N+2) group 

We consider a fermion system with N single-particle states. Let c a and c£, a = 1, . . . , N, be 
the annihilation and creation operators of the fermion, respectively. The SO(2N+l) canonical 
transformation U(G) is generated by the fermion SO(2N+l) Lie operators consisting of the 
set {c Q ,4,£« = c t C/3 - i.S^Eap = c a cp,E a ? = etc*; a, (3 = 1, . . . , N}. The U(G) 
induces an inhomogeneous linear transformation (TR) (for zj^l) in the space spanned by 
the fermion annihilation and creation operators (not a linear one) involving a g-number 
gauge and is specified by an SO{2N + 1) matrix G as 



U(G)(c, c\ -^)U ] {G) = (c, c\ ±=)(z-p)G 



G = 



A B* -X* 
B A* X 
Y -Y* Z 



X- X Y- V 7 



G^G = GG ] = l 



2N+1, 



U{G)U{G') = U{GG') , U{G- l ) = U- 1 {G) = U\G) , £/(!) = ! 



(2.1) 

(2.2) 
(2.3) 



where (c, ct,4=) is a (2A^ + l)-dimensional row vector [(c a ), (c^) , 4=) and A = (A") and 
B — (B ai ) are NxN matrices. The symbol * denotes the complex conjugation. The operator 
p in the gauge factor z—p is defined as p = x a c { a —x k a c OL and satisfies p 2 = — x^x a = z 2 — 1. 
When .2 = 1, then G becomes essentially an SO(2N) matrix g given by 



a b* 
b a* 



9*9=99* : 



'■2N- 



(2.4) 



The HB (SO(2N)) WF \g> is generated as \g> = U(g) |0>, where |0> is the vacuum 
satisfying c a \ 0> = 0. The WF | g> is expressed as 



(7> = <0|t/"(^)|0>exp(J-g aj8 44)l > > ? = &a 1 



<0|C/(p) |0> = [det(l-g*g 



4 2 , r = - • In 
2 



det(a* 
det(a) 



(2.5) 
(2.6) 



where det represents determinant and the symbol t denotes the transposition. By a = (a") 
and b= (b a i) we denote NxN matrices satisfying the orthonormalization condition for the 
HB amplitudes a and b jl] . By q a N x N matrix is denoted which is a variable of the 
so ( 2Ar ) coset space. On the other hand, the SO(2N + 1) WF \G> is generated as [23] - 



U(N) 



G> = U (G) | > and can be expressed as 



G> = <0\U(G)\0>(l+r a cl)exp(^-q al3 cl e cl)\0>, r a 



1+z 



y X a ~\~qa/3X j g) 



<0\U{G)\0>- 



1+z 



[det(l-g*g 



4 e l 2 



(2.7) 
(2.8) 
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The SO(2N+l) group is embedded into an SO{2N +2) group. The embedding leads to 
a unified formulation of the SO(2N+l) regular representation in which paired and unpaired 
modes are treated in an equal manner [23]- [21]. We define (A/+l)x(A/+l) matrices A and B as 



A= 


A 

y 


x* ' 


, B = 


' B 

y 


x - 




.2 


2 . 




. 2 


2 J 



A = a- 



x y 



2(1+ z 



; B = b + 



xy 



2 1 + z 



y = x T a-x ] b. (2.9) 



Using the orthonormalization H2.2[) . the matrices A and B are shown to satisfy the orthonor- 
malization condition for the (jV+l)-dimensional HB amplitudes and to form an SO(2N+2) 



matrix Q. The representations for Q and for ^f^p^ coset variables Q are given as 



A B* 
B A* 



SO(2N+2) 
U(N+1) 

, &g = g& = i 2N+2 , q=ba~ 1 



q r 
-r T 



(2.10) 



which shows that the SO(2N+ 1) variables q a p and r a are just the independent variables 
of the ^fn^-f) coset space. The paired mode and the unpaired mode variables g Q( g and r a , 
respectively, in the SO(2N+l) algebra are unified as the paired variables in the SO (2N+2) 
algebra. We denote the (iV + l)-th dimension of the matrices A and B by the index and 
use the indices p,q, . . . running over and the single-particle states a, (3, . . .. 

Expectation values of the fermion SO(2N+l) Lie operators, i.e., the generators of rotation 
in (2iV+l)-dimensional Euclidian space, with respect to \G> are given as 



<£ ,Q g + -(5 Q/ 3 >G — Raf3 



<E a p >g 



l -(B* m B Pl -A« l A^) + l -5 a ^ 
K ap = i (AFfBpt - B m A^ , <E^> G = K 



<c a > G = K a0 -R i 



ao-^i^tyi+B^i 



<ci> G = K* -R, 



a • 



(2.11) 



The unified matrices Rp q and K pq are expressed in terms of the ^f§j|pp coset variable Q pg as 

R P9 = -[Q k Q{l-Q k Q)- 1 ] vq , K pq =[Q(l-Q*Q)- 1 ] pq . 
The expectation value of a two-body operator is given as 



pq 
(2.12) 



Let the Hamiltonian of the fermion system under consideration be 



H = h afj [E a p +-S aP \+-[aP\jS)E^E, 



(2.13) 
(2.14) 



which is expressed in terms of the generators of rotation. We call such a sytem a fermion 
top. The matrix h a p related to a single-particle hamiltonian includes a chemical potential 
and [a/9|7(5] = — [a<5|7/3] = [7^ a/3] = [/3a|$7]* are anti-symmetrized matrix elements of an 
interaction potential. Parallel to calculations by the usual HB factorization method (See 
Refs. [2] and [3]), the expectation value of H with respect to \G> is calculated as 

1 

— < 
2 



<H> G = hapKE 01 ^ + -5 a p >g 



+ -[aP\ 1 5]l<E 



-U a p > G < E\ + U^ 8 > G +-<E^ > G < E S p > G 



(2.15) 



2 M 2 " 2 

Then the <H>g (I2.15P represents approximately the energy of classical motion of the 
fermion top. 
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3 Boson images of fermion SO(2N+l) Lie operators 
and equations of motions for bosons 



According to Fukutome [23]-[23], the fermion SO(2iV-|-2) Lie operators {E p ' E pq , E pq } are 
constructed by the fermion S0(2N) Lie operators {E a p, E a p, E a/3 } and the operators 



E 



E< 



aO . 



C 1 " P,= 

±P E° 



-P-cl P± = -(l±(-l)»), 
= c a P+, E° = -(P_ — P + ) , E a0 - 



aP- =P+c n . 



(3.1) 



The boson images e p etc. of them are given in the following forms pH], [TT] : 



d 



pr 



OB 



-B, 



d , ^ d d 

. Al* u A p 



qr 



d 



£pq AP *dB 



■B, 



qr dB pr 
d 



r 8A P r 



■A 



q* 



r dA q r 

d 



^-+B — 

dBL pr dA\ 



' pi 



' <3' 



£ t =-£* 

pq pq J 



(3.2) 



qr PiA p r FiK* P r a ai " ' "pq "pq' "pi "9V 

■> q r yJJ~\ r pr 

Using ( 13. 21) . the images of the fermion S0(2N+1) Lie operators are expressed [23]-[24] as 



which are rewritten as 
1 



■■£ a)3 C 



-=-0q t a , c a — fc fc 0' 



(3.3) 



^[C a ,Cp\ 



A a jBp ? - AVb^, E? p = -E. 



a0~ 



JTtk 

'■^api 



■■y/2(A%yr+y\pJ) ,yr=^=( 



V2 

X a i -^N+l —^on B a Q = X a , B 



*A- B n ~ 



aN+l ■ 



yr+N^-y*, y Q =z, y N+1 = -z,(r=o, 



Or J i a' 

N,r=0, 



'•t — — /-** 



> (3-4) 



N,N+1,...,2N+1),) 



where A p q and A p q1 etc. are defined by all the variables A p q and A p q , etc., in Q (I2.10p and 
their partial differentials as 








q dAZ 




A* P T 







[A^,A^ s ] = [A^,A T \] = 6 pr 6 qs , A*> q \ 



-A 



*p 



9 - Q A P I ' 



= 0, 



(3.5) 



[A\, A\] = [A p q , A Tr s ] = 0, [A* , A Tr J = [A* , AM =0. 

Similar definitions hold for B in order to define the boson operators B pq and B* q , etc. Here 
we have used the same symbols as those used in Refs. [10] and [TT] . 

The SO (2 N+l) Lie operators are the generators of rotation in the (2iV+l)-dimensional 
Euclidian space. The fermion top has two kinds of coordinate TRs, namely, the TR for space 
fixed coordinate frame (left TR) and the TR for body fixed coordinate frame (right TR). 
To discuss the c-number limit of the bosonized fermion SO(2N+2) Lie operators, we make a 
coordinate TR bringing the body fixed coordinate to the frame referred to a quasi-particle: 



A p ~=A p ~0^, B pq = Bp 7 0^ (q,r,= 0,...,N,N+l,...,2N+l), 

O r p = A p, O r+ N+i,p = B r p } (boson operators in the quasi-particle) 

p * p ~* ~p ~ ~ 

r+N+l — /J pr' — ^pi "^-N+l—^pi °p,r+N+l 

Y r+ N + i = —y r , y =z, y^ + i = —z. 
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p* 

r ) O p0 : 



-X^B 



Xpi °p,JW : 



W3.6) 



Though we use the same symbols parameters, A p ~ etc. in (13. 6p are not identical with those 
involved in the variable Q. They are coefficients of a time dependent coordinate transforma- 
tion for the Q. The coordinate transformation is represented in the matrix form as 



(3.7) 



The bosonized quasiparticle SO(2N+l) Lie operators with indices i and j (i, j = 1, . . . , N) 
are constructed from the operators in (13.61) in the same way as (13. 4p and the operator W~ 
is introduced in the following form: 

l r , t - sits 



r ap 

q 


B* ' 

pq 




Ap 


B* 

^pr 




SI q 


B r q 


B 

. "pq 


A p * 
q - 




B 


AP* 




B 


^ q . 



E^-idl^B^-A-Ar, 



Eij ^ [di, dj\ A. 7f. Bj r JK. ~ Bj r , E 1 ^ 



4 K- 



d l =V2{ 'A l \yr+ylB l r) , y F =l(4-»Sr), 4 





' B ' 






A p 


. > 




y p 





} (3-8) 



In the boson images of SO(2N+l) Lie operators, particularly the images c a and ct must 
satifsy the anticommutation relations {ct , c^} = 5 a/ 3, etc., when they operate on the spinor 
subspace. Therefore, we use in the regular representation space the following image of the 
Hamiltonian with the Lagrange multiplier terms to select out the spinor subspace: 



H' = ht aP {{cl, c }-5 a p)+-l a p{cl t , 4} + i^{c Q , c^}, (k* a p = kp a , l a p=lp a ), 



(3.9) 



Using ( 13. 9p . ( 13. 5 p and ( 13. 6p . the Heisenberg equations of motions for W| (p = 0, . . . , N, N+ 
1, . . . , 2N+1) are calculated as 



^{F aP ,B^}+^{D aP ,A%}+^={M a ,y ¥ } 
-l{Fi,,A%}- 1 -{D^,B,p} + ^{Ml,yp} 
-^={M a ,A a p} + ^={M^ a ,B a p} 



(F a pBpq+D aP A $ ~+ V2M a y^j Ogp+-{f ap B P r+ d aP A%+ ^2m a y~ r , C%} 
{-F* p A%- D* a pB$q+ V2M*y«) Ofp+ \ { -fJpA%- d^Bpr+ ^tn^, <%} 
(V2M a A a ~+V2M*B a g) 0^+-{V2m a A a ¥ +V2wiB a r, <%} 



(3.10) 



where the W~ denotes a classical part of the W~ and the F a p, etc. are defined as follows: 
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F^hap+WM (E^ + U^F^+f^ f a ^[a^8] {E\-<El> G )=f^ 



X3.ll) 



M Q = k aB c B +l. 



The SCF parameters F = (F al3 ) = F\D= (D a p 



Kb {c-b -<cb>g)+Lb (4 ~ <C /? >G ) • 

-D T and M= (M a ) are defined by 



7 aB = h a/3 + [al3\>y5}R~ f s, D a p = -[aj\(35) 



K Sl ),M a . 



'■k a B<CB>G + l a B<c\>G ■ ( 3 -l 2 ) 



which involve effects of unpaired-mode amplitudes as well as those of paired mode. 

The Hamiltonian (I3.9P is expressed in terms of the boson operators, i.e., an assembly 
of oscillators, but represents a (2iV + l)-dimensional fermion top. From (I3.10p the classical 
equation of motion for the W~ of the fermion top is obtained as 

V2M a 
V2M* 



iHWS=T a pW§, F a 



F a B 

-D 



D 



aB 



aB 



aB 



y/2Ml V2MJ 



(3.13) 



Simultaneously the direction of the axis of the rotation fluctuates quantum mechanically 
around the classical axis. This picture of the fermion dynamics was first given by Fukutome, 
Yamamura and one of the present author (S.N.) [6]. Equation ( 13 . 13[) is transfomed to 



i^W a = F a0 Wf and ih 







- Xq - 


71 




A 


x a 


— FaB 


X B 


V2 






z 




Z 



(3.14) 



using the relations in (12. 9p . from which we obtain the equations of motion for HB amplitudes 
a and b though including unpaired mode amplitudes x and x*. 



iha- 
ihb- 



l + z 



(x*M T + M*x r ) 



a — 



D* - 



l+z 



(x*M j! -M*x f ) 



■(xM t + Mx t ) 



b+ 



D 



l+z 



(xM T -Mx r ) 



(3.15) 



a. 



This is an extended TDHB (ETDHB) equation in which the paired and the unpaired modes 
are treated in an equal manner [6] and is applicable to both even and odd fermion systems. 
Then it is the extended TDSCF theory having the same level of a mean field approximation 
as the usual TDHB. By putting a and b as e~a and e~b, we obtain a static EHB equation 



1 



l + z 
1 



(x*M T + M*x T ) 



1 



l+z 
1 



(x*Mt-MV) 





a 


= E 


a 


,F 


- x - 

A 

X 




b_ 




b_ 




"71 

z 



-0. (3.16) 



Though we use the same symbols, a, b, x and z in ( 13 . 16[) are time independent amplitudes. 
The constraint term in (13. 9 p should vanish in the physical fermion space but their classical 
part might not. Using the third of (13.121) we must determine the parameters k and / in M. 
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4 A Solution of the extended Hartee-Bogoliubov eigen- 
value equation 

We use a spherical symmetric single-particle state specified by the set of quantum numbers 
{n a ,l a , j a ,m a } which is denoted as a. The time-reversed single-particle state a is obtained 
from a by changing the sign of m a . We use a phase factor s a = (— l)-? a_ma in the time reversed 
quantity. The SO(2N +1) density matrix has the time reversal symmetry R^p = s a spR-p. 

By introducing a new amplitude b defined by h m - = s a bai [28] and using the time reversal 
symmetry F* j3 = s a spF^p, the EHB eigenvalue equation is transformed into 

f ap-J^ z (<Mp+M*xpyfa ^{spDlp- ±- (xlspMj-MZspxtybp^Eta* ^ ^ 

(4.2) 



(4.3) 



_ _ (s a XaSpMj+s a MaSpX^bp—^s a D & p - (s a XaMp-s a MaX p)^ a 13 — E { b, 

~ F ^("^l)"^(^l) + ^ M ^ =0 ' ^( s ^)+ s ^(-^|)+v / 2^M^ = 0. 
The SCF parameters F a p and D a p are decomposed as 

F a/3 = f a p + ^j^[al3\-fS} (x*<c s > G +<ct f > G x s )=f a p+-^, f j x = f x , 

Da f ) = da P + 2{l+z) \ [ a7 '^ ( x j <c s>G-<c 1 >GXs)=d a p + -^ 1 d^ = -d x , 
where f a p and d a p are the usual SO(2N) SCF parameters defined as 

f a p = (e a -X)-5 a p+[al3\-i5}p 1 5, (/* = /), d a fi=-[aj\fiS\(-K Sl ), (d T = -d) } (4.4) 
and the SO(2N) density matrices p a p and K a p are given by 

p a p = <E a p + -5 a p> gj K a p = -<E a p> g , K* a p = <E a P> g . (4.5) 

The expectation values of the annihilation and creation operators, with respect to the state 
\G>, are given in terms of the matrix entries p a p and n a p and the amplitudes x a and x* a as 

<C a > G = ^S a p-p* a p^jxp + K a pXp, <C j a > G = <C a >G. (4.6) 

We solve the EHB eigenvalue equation by the two-step diagonalization method [12]. First 
we diagonalize the separable particle-hole-type interaction and neglect its exchange effect. 
Then, the Hartree approximation leads to the eigenvalue equation 

t a pWpi = eiW ah WaiWaj=5ij, W* ai Wpi = S a p. (4.7) 

The time-reversed Hartree state t~-w a j is degenerate with the Hartree state w ai in energy: 

fU t ~i w P~i =e ~i t ~i w oTi^ e ~i= e i- (4-8) 
The phase £j, the counterpart of s a , satisfies the relation t~ { = — The Hartree state w ai 
satisfies the relation w'^ i = s a t-rw a --. We put the HB amplitudes into product forms 

a"=Mjiy Q i, b a i = ViW a i, (not summed over i) , \ui\ 2 + \vi\ 2 = 1, (4.9) 
for which we impose also the time reversal ansatz 

s a tia a -r = a* a i7 u 1 = u* J s a ti6 a 7 = 6^, vr=v*. (4.10) 
We also introduce the following quantities associated with the unpaired mode amplitudes 

Mi = M a w ai , Mi = s a M^w ai , Xi = x a w ai , Xi = s a x* a w ai , (4.11) 

which have the time reversal properties t%M~ — Mi, t%Mj — Mi, t%x% = Xi, t~x% = Xi. 
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Further assume the pairing potentials d = (d a p) and d x = (d x>a p) to be constant. This 

makes the situation very simple as in the BCS theory does [29]. Then, the pairing potentials 

have the following forms: 

1 1 

d a p=s a 5 a pA, A=-gs 7 K^ and d Xj0 ,p=s a 5 a pA x , A x =-gs 7 (x 7 <c^> G -<c^> G Xy) , (4.12) 

where g is the strength parameter for the pairing force. By substituting (14. 9p into (14. ip 
and using the Hartree equation (14.71) and the definition fx,i = w aifx,a/3 w i3i( = fx,i)^ the EHB 
eigenvalue equation ( 14. ip under the assumption ( 14. 12ft is converted into the following form: 

eiUi-^r— (x*M i +M*x i -f Xii )u i +A k v i - TTZ (i-¥,- — M*Xi — A*) v i = E i u ii 



1 + z v , . , . . . 1 + z 

-e i v i + T — [x+Mi + M^Xi-fxA Vi + Am + 



1 + z 



1+z 



x*M i -M,x i + A x )u i = E i v i . 



(4.13) 



Here we have used the orthonormalization condition in (14. 7p and the new quantities in (14. lip . 
Similarly, multiplying ( 14. 2 j) by and summing over the index a, we obtain the equations 



1 

72 



eA 



fx 



1+Z, 



x, A4 



A 



V2\ 1+z, 



i+z, 



x 



rWlzM^O, (4.14) 



from which we have the formal solutions for Xj and X; 



Xi 



Ef 



2z 
E? 



eA 



fx 



1 + Z 
2 



MA- A4 



A, \- 



1+z 



■)Mi},Xi = - 



2z 
E? 



eA 



fxA 



1+Z 



Mi- A*-+ 



A! 



M 



1 + z 



+ 



. A x 
A+ 



1+z 



► (4.15) 



In the second step we diagonalize the pairing interaction contribution only to the SCF 
pairing potential D = (D a p). In order to diagonalize the EHB eigenvalue equation (I4.13p . it 
is very useful to notice some important relations satisfied with x^ x^, Mi and M». Using the 
above formal solutions (14. 15ft . we obtain the relations 



2 z 



fx, 



x*MAM*Xi=- 1| <i 2 ( e 

_ 2z 
x*M-M*Xi=-^- 

E 2 



x *M+M* Xt = h k+gM \M\ 2 +{ A+^-j M*Ma-( A'+^J 



A„\ . .._ 



AtW 



fx 



1+z 
A* 



Mil 



A+^ 

1+z 



M*M~ A*4 



2^ 



A* 

"l+Z 



> (4.16) 



A*+^J(|M| 2 +|M i | 2 ),^M-M J x l =-^^A+^J(|M| 2 +|M| 2 ) 



A.., 



Substituting (14.161) into (14.131) . we have a secular equation 

D 



Dm 



'm 

— p-—F^— E 

* M t 



(4.17) 



where the quantities Fu and D M are defined as 




M*Mi+ ( A*+^J M, }> + f x 
M*Mi-[A*+^)M*Mi}+f x 



> (4.18) 



From ( I4.17P and f 14 . 1 8 j) we obtain an eigenvalue 
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Ei=Ei + 



2z 



2 +lMI 



+ 



1z 



fx 



1+2, 



l+z 1+2 

which reduces to the usual HB quasiparticle energy if z = 1 [2]. It involves unpaired-mode 
effects in contrast to the usual HB theory, which is unable to describe the unpaired modes. 
With the help of the EHB eigenvalue equation (14.131) . the eigenvectors are obtained as 



fx.i . 2z 1 f / f x 

2 e,-' 



U; 



l E&Gi] l+z ' l+ Z E 2 



1+2 



A*4 



a: 



1+2 



.1 /,.!/, 



2z 1 



e ? 4 



/a; 



1+2 



(|M| 



|M 4 | 2 )+( A+^jJ^M I +[A^ 



Ei-e 4 , 



fx 



— z^M- ) |M| 2 +( A+^-]m*mM A*+^W 
l+z l + zE 2 \\ l+z)^ 1 V l+z) h \ 1+2 ' ' 



> (4.20) 



2^ 1 



e 4 +^)(|M| 2 -|M| 2 )+(A4 



1+2/ V 1+2 



MM 



1 + 2 IV 1+- 2 

From (I4.15p . we can get the solutions for the unpaired-mode amplitudes and for z as 

1 1 n f e _|_ fx,i 1 i « r H i i a , 1 9 ilTT 

If 



3 i 



l+z 



IM 



IA+ 



1+2 1 



IM, 



+ e 



1+2 



A+— ]M*Mi+(e+ ^" 



1-12 4 ^Jf^ 
Et IV 1+2 



1+2 



IMJ2 + IA+ 



1+2 



A*+ 



1+2 



A, 



1+2 ' 



e+ 



1+2 
1 



A4^ 



MM"-ie 

1+2/ \ 1+2 

1 



IM 



A*+^JM*M 



1+2E^(IM| 2 + |M,| 2 ) 



(4.21) 



E? 

where we have used the normalization condition J7 |xj| 2 = ^ i |xj| 2 = l — z 2 . 

Up to the present stage, the SCF Hartree potential, the chemical potential A and the 
pairing potential remain undetermined. With the use of (I4.3p ^( l4~6l) and the product form 
of the HB amplitude (14. 9ft . the SCF Hartree potential given through (14.31) is determined as 

fa/3=(£a-X)-5al3j2i \v i \ 2 W ai W^ i +[af3\j5] J2i \Vi\ 2 WyiW^ 1 

if 1 1 1 X4-22) 

fx,ap=2 \ a P 1 7<*] E ify ( 1-1 Vi 1 2 -| vj | 2 ) x*Xj+- (upj+UjV*) x*xj+- (u+Vi+UivDx+Xj I w 7i w* 5j , I 

and the pairing potentials and the chemical potential are determined by the conditions 

& = lgEi(utvi+ Ui vt), & x = -\gY Ji {<v i + vtu i )(\x i \ 2 + \x l \ 2 ), and (4.23) 



■Af—Raa—J2i 



i I Vi | 2 | Xi | 2 - j (u*Vi + UiV* ) (x*Xi + X*Xi ) 



(4.24) 



where M is a number of fermions in the system. They contain unpaired-mode effects and 
reduce to the equations for the gap energy and the chemical potential of the usual HB theory 
if we put 2=1 [2]. The magnitude of the gap energy evidently decreases due to the presence 
of the unpaired-mode effects (Pauli blocking effects), which is also seen in the solution of 
the improper-HB eigenvalue equation presented by one of the present authors (S.N.) et al. 
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5 Energies of classical motion and quantum mechanical 
fluctuation of fermion SO(2N+l) top 

Through the coordinate transformations (13 .6p and (13. 7p . the fermion SO (2 N+l) Lie oper- 
ators are expressed in terms of their quasiparticle expectation values (c-number) and of the 
quasiparticle SO(2N+l) Lie operators (quantum mechanical fluctuations) as follows: 

ETp^E^^+W^E^^+U^-BgiBX^E'^i^+D^A '*E lj +A" l D r E 



'13 



+ \ faxp + d l+ l - <A%x p + d\ , 



E 



a/3 



< E a/ 3 > G +E a /3 = < E a p > G +[B ai A^j—Bi3iA a *)[E l j +-5ij}+A a *A p *Eij+B ai Bi3jE 



+ ^ (A a txp - A^x^j di + 7; (BodXp ~ B^x a ) d\ 



' i j 



K5.1) 



c a = <c a > G +c a = <c a > G - (B ai y*+ A a *y^ [E )+-5 l3 ) -A a *y*E lJ +B ai y J E 



+ ^A a *+-x a !,:)d l -\ [zB„ 



^x a yi)d\. 



The combination E l j+-Sij, in which the constant term — -8^ in E l - is eliminated, is relevant 

to describe the quantum mechanical fluctuations. 

Using (15. ip . the boson image of the Hamiltonian with the Lagrange multipliers (13. 9p is 
transformed to the following form: 

H=E + H, 



E o = <H> G + k a/3 [<cl>G<Ci3>G~S a ^+^l a ^<cl> G <cl> G +^ <C ct > G <Cf3> G , 



H 

1 



1 ~a 1 ~af 1 ~o.fi 

■^F af} E p + -F* p E p + -D a pE - 



1 



-D*oE a p+M a c r a +M*c 



>(5-2) 



uP\l8\[{E a,E 5 }+-{E , E sp } +k a $? a , c p }+k*g {H , c a }+l a(3 , cU+Z* {c a , c p ) 



where Eq is the energy of the classical motion of the top and deviates from the quasiparticle 
expectation value of the Hamiltonian < H >g owing to the presence of the constraint term H' 
in (13. 9p . The coordinate transformation method contrasts with the coherent-state method 
[32] . The Hamiltonian H is the boson image of Hamiltonian for the quantum mechanical 
fluctuation. To treat such a quantum mechanical fluctuation, we make the successive coor- 
dinate transformation Q similar to (13. 7p . We may decompose the tranformation matrix Q 
approximately into a product form of two matrices in the following form: 



B 



~ ~Tp-k 



a* a 



pr 



Bp r A P r 



A r „ B 



B rq 



rq 



Ar-k 
^ q 



A p B* 

q ^pq 

B A p * 



A p B* 



Bp r A p r 



2l r „ 03 



rq 



q . 



, (5.3) 



in the second equation of which, the matrix (3 stands for the fluctuation matrix deviated 
from the statioary matrix Q . The decomposition of the matrix Q (15. 3p first proposed in Ref. 
[31] is equivalent to making the Bogoliubov transformation by the stationary amplitudes A 
and B that brings state vectors and operators into the quasiparticle frame determined by 
the stationary solution. From the decomposition we have approximate relations between the 
amplitude in the Q quasiparticle frame and the amplitude in the (5 quasiparticle frame as 

J^P q J\P ,p 5i g~\~13pY,^&rq cllicl 13 pq 13 p7-3i H-^^ 'y-^rfj* 
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The bosonized quasiparticle SO (2N+1) Lie operators in the Q quasiparticle frame with 
indices % and j (i, j = 1, . . . , N) are also constructed from operators similar to the operators 
in ( 13.61) along the same way as (13.81) in ^the following forms: 



E, 



2 1 
1, 



'■'./ 2 A r 13 jr A ~ 13i r , E E^j -Ejjj 



V tl 



(5.4) 



Using (I5.3p . the fermion SO(2N+l) Lie operators describing the quantum mechanical fluc- 
tuation are expressed in terms of the quasiparticle expectation values (c-number) of them in 
the Q quasiparticle frame and the new quasiparticle SO(2N+l) Lie operators as follows: 



E 



< E% >r3 +E%=< E% >^+{A%A^-B m B*J [E\+-5^ m A%,+A^B 07 E 



+~[B*M+A p :x* a ) d^~ [A^+B pi xl)d\, 



E a f3 =< E a8 > G +E a p =< E aB > G +[B ai A^—B Bi A a *\ \E l j+-Sij \+A a *A p *Eij+B a SpjE 



; ij 



+ ~[A a *xp-A p *x a ) d i + : - [B ai x B -B Bi x a ) d\, 



W5.5) 



c a =<c a > 5 +c a =< Ca > 5 - [B^yp-AJyi] [E —A a *y*Eij+B ai yjE « 



+ [zA a *+-x a y* di+[ zB ai --x a yi d 



A~,t 



where the quasiparticle expectation values in the Q quasiparticle frame are given as 



I . ~^ _ 1 



<E a + -5 aB>d =R a 



<E aB >g 



-K n 



B*Ba i ,-A"A p :)+-6 l 



- (A a *,B S i' — B a i/A^,) , < E a/3 = — < E aB >% 



■>/3i> 



'at' ■"- i' \ i 



G 



G' 



<c a >c= ^(A a pyi>+B ai/ y*) , <^ >Q =<c a >*~ . 



G 



(5.6) 



whose forms are the same as those of (12. lip . Substituting (15. 5 j) and (15. 6 j) into (15. 2p . the 
fluctuating Hamiltonian H with the Lagrange multipliers is converted into 

H=E~o + H 



r 



l 



l. 



[a0| 7 5]({E% EP^E-\ E sp }}^k a p{cl c^k*^, c^l* a0 {c a , c } 



E Q =)f a p<E^>Q^F*e<E*>a^D a p<E^^ 



+- A [a^5] [2 <E^> d <E\> d + <E*r> d <Esp> d 

+ k a p<cl> d <Cf3> d +-l a /3<cl>G<C i p>G +2^ <5q> G<C/3>G: 

+\ [aP\j6] (<E\> 5 E" + <E" 0>G lF s +\<E sp > G H^> 6 E sp 



X5.7) 



+-k aB (<c b >g 4+ <ct>^ c fl ) +h.c. + T Lb (<eti>;5 4+ <c£ >« c£) +h.c 



*G c /3 



^oc G 



4 -a/3 ^/3^G c «" 

where h.c. denotes the Hermtian conjugate. A part of the fluctuating Hamiltonian H is 
given up to the first order in the boson images E and c. 
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The part of the fluctuating Hamiltonian iifis devided into two parts as H= H~.+H^. Above 
~_ ~_ d hi 

all the is important. H-^ is given as 



H 



+ 



eL+h.c. 



dj+h.c. 



+i [a/3| 7 5] 



i 7 



eL+h.c. 



16' 



[aP\j6} <Esp> d (B^-mx* a ) - <^> g (A s *x/3-A^*xs 



dj—h.c. 



(5i 



+ 



d, 



^+h.c. 



+ 



<4>G MV^a2/J+ <^>G 



d^+h.c., 



2 aai.i 1 "a: (j ( 1 2 

which is exactly the so-called dangerous term and should vanish. As for the annihilation 
opertor dj/ of the quasiparticle state i', the right-hand- side of (15. 8 1) is calculated as follows: 



F Q/3 ^x^fi^j+F^ [A^xfrB^xj-D,. ^B^-Bl^-D: ; [A<;,r ; -A ;,.?•„ 
^M Q 4M* (zA«*+±x a y*, 



+2k a p 



< E %>6 [B^xs+A^x^ ~{<E sp > d (B* tf x;-B; tf x* a )-<E^> d [A%x p -A%x & 

t 1 



= 2 
+2 



L 



<4>G (^ai'-7?ayt^<$a>G W^'lf^ 



c^a^vy a \ 2 



x*+2 



+1, 



<C I3>G (zA a -^ a yt^<Ca>G W^i^Vv 



f* b%,+d^ap;,+V2m*(% 



V2J\ 



V2M a B* m , + y/2M*A°^\ V2z 



(5.9) 



2£«, L4^x*-^,x a -5^J +V2F (y^M Q E* 4 ,+ x/2M*A 



+ 



FapA^+DapB^, + V2M a f- ^= 



x* + 



?7? 
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The vanishes if we use the normalization condition A a px^,— B^ ir x a — zy*, = which is 
derived from ( 12. 2p . We define the tilded SCF parameters in the G quasiparticle frame as 

F a p = h aj3 +[aP\"/6]R 7 s, D a p = -[aj\P5)(-K Sj ), M a = k aP <Z p > 5 +/ Q/3 <cj> g . (5.10) 

In (15.91) we also have used another eigenvalue Ey related to a quantum fluctuation given by 
the same type of the eigenvalue equation as (13.161) . Namely, the eigenvalue equation for the 
fluctuating Hamiltonian with the eigenvalue E# is given as 



p 



L v^J 



■Ei> 



L V2\ 



a/3 



V2 



x 



p 



V2 



-0, T t 



a/3 : 



a/3 



D 



-D 



a/3 



afi 



a/3 



V2M a 



V2Ml V2MJ 



(5.11) 



The Hamiltonian .H— describes the quantum mechanical fluctuation and is given by 
hi 



H 



E 

A a lF a ^;+D a ^*+V2M Q 



y 



bJf^F.+d^b^m* 



^)\-B ai [F^ J+ D^;+V2M: ^ 

y 



A^FapAt+DcpBfc+y/ZMa (- 



^+AJ^B; i+ D*^:W2M- 



E 



'.i 



E ij 



-y^M^+M+A 
1 

+ 4 



la* 

3 



(5.12) 



B* , <| F a ^; + D aP B* pj + V2M D 



yl ( M a B* ai + M* A° 



+Aj\F^ i+ Dl p A^W2M^ 



E 



B^F^+D^+V^M* f--^)} +A^F aP Bp j +D aP A p j +y/2M a (j=) 



+y 3 [M a A%+M*B c 



E ij . 



Using eigenvalue equations (I3.16P and (15. lip with eigenvalues E^ and Ei and normalization 
condition G^G — GG* = la/v+ij the Hamiltonian LTjj, (15. 12ft is simply rewritten as follows: 



(5.13) 



E - E j [A«AJ + B ai B* aj ) (E )+-5 t j +-E 3 M^A? + B ai B* aj + -y^j [E 3 +-5 l3 
+ 1 -E J (B^AJ+A^B*^ E l3 - l -E 3 (B m A° J+ A«B^Ev 
+- A E j (b^AJ + A^B^ - Ijffig) E l3 - -E 3 fe»2« +A^B a] - l -y t y 3 ) E 



, 'J 



^[Er^EA IE \ —E^yAE )+^5 l3 \ —E&yAE ^M^tfffiE^tfjftS «. 
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Now we consider an excitation energy due to Hj^ (= H) (I5.13P in which y and y* are relevant 
variables to treat the excitation energy. This contrasts with the role played by the variables 
x and x* in (I3.16p . We attempt to derive the excitation energy by the Sawada's eigenmode 

method [13]-[T3]. The method, however, does not work well because the Hamiltonian H-^, 
contains the unpaired-mode amplitudes. We return to the traditional Bogoliubov's proce- 
dure [3]. Let us introduce the following annihilation and creation operators di and d\: 



di, d\ 



d u d\ 



Ui Vi 



-Vi Ui 



di, d] 



9i, 9i9i = 9i9i = h- 



We also define a new bosonized quasi-particle SO(2N) Lie operators as follows: 



(5.14) 



(5.15) 



Ey = -\d\,d 3 \, E\ = E% E^-^dj], E' 1 -E u E' r 

Through the transformation (15.141) . the old bosonized quasiparticle SO(2N) Lie operators 
(15 .4p are converted into 



1 



E) = ^[-v*d l + Uid\, u* dj + vj d) ] = uiu* E)+v* Vj E 1 }- v* u* E^ + urfj E*, E\ = E\, 
E ^ = Ztf* E ) +^,*u j E*j- u?S* Eij +%^j E * , 



<5.16) 



E ij = ^[ u * d i+ v id\, ") d i ■ r j d j. =ViU*E)-u*VjE t ]+u*u*Eij+ViVjE 



it 



E^ = -UiV^Ei + v^UjE^ + v^Eij + UiUjE^. 

Substituting (I5.16P into each term in the last line of (15.131) . we obtain the Hamiltonian H-jj, 
in terms of the new bosonized quasiparticle SO(2N) Lie operators as follows: 



— A Ejy*yi[ UiU* E )+v*v j E i j-v*u*E ij +u i v j E» + -5.. 



4 
1 

T 



— E { y\yj ( ViV* E ) + u* u 3 E'j-u* v* E tj + v { Uj E^ + ^ 5, 



E J y*yt(v i u*E)-u*VjE i ; +u*u*E ij +v l v j E 



'.i 



— E jyj yA - Ui v* E V + v* Uj E'j+v* v* E vj + u iUj E 



(5.17) 



1\ 1 



- \ ( E jyjy^t^ E ivtmt^ (e '] ■ ^ 

+\ - \ (^jViUi-Eiytv^ (^jUj-y*v^E^. 

The last two terms correspond to the so-called dangerous terms and should vanish. Then 
we can get the relation between u* and v* and consequently have a simple and important 
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result for the state j and also the same for the state i 



u *i = % v *v (";".' • r ) r :> 1 



UjUj = v j v j = -, u i u i = v i v i = -. 



(5.1* 



Substituting (15.181) into the second and third lines from the bottom of (15.171) . and putting 
Vj = Vie t9 which is due to ( 15. 18ft . then we have the Hamiltonian Hjj, in the following form: 



1\ 1 



2/ 8 



id I -p i] _ 



The Hamiltonian i2~ can be diagonalized as Hj n =^2 i Ey i \E\+-} where E^ is a diagonal- 



E 



ized energy which is given, e.g., for i, j = l~2 as 
1 



Er, 



E 1 +-E 1 )+ (E 2 +-E 2 ) + J{(E 1 +-E 1 )- (E 2 +-E : 



■(Ei-E 2 )\yi\\y2\sm9 



(5.20) 



The diagonalized energy Eyi is given also in the case, i,j=l~N but then it becomes complex. 

Using (15. 2p . (15. 7p and ( 15 . 1 9j) . finally we obtain the boson image of the Hamiltonian H 
with the Lagrange multipliers in terms of E (energy of classical motion of the top), E 

(energy of fluctuating motion) and quadratic terms of fluctuating operators E and c as 



H=E + H 



1\ 1 



Eo + Y,ME\ + - + 7 M 7 i] [{E%El}+-{E^,E tp } 



2 4 



1 1 11 



>(5.21) 



, 1 



= E + E + j: i Ey l [E\+-j + -[ a /3\iS\U& p ,& s } + -{E»r,Esp} 

£ = <#>g+^ (Afa /8 <4> G +M*<c a > G )-^A; aa , 

+ h a / 3 <E a l3 > G +h* a!3 <E a ^> G + (F aa — h aa )^ , 
E ^(F a e<E«p> s W^<E a t> d +D a ^ 



F a p<E a „> d +F*s<E a>c+D a p<E a P>Q-D* a p<E a p>^ a p<E a p>^* oi0 <E „>g 



>(5.22) 



+ -{M a <& a > s +MZ<c a > s y r 

The components of the generalized angular momentum, i.e., the fermion SO(2N +1) Lie 
operators, are not in general commutable. Therefore, owing to the presence of the fluctuating 
Hamiltonian H, the direction of the axis of rotation fluctuates quantum mechanically around 
the classical axis determined by (I3.13p . If we consider the quadratic terms of the fluctuating 
operators in (I5.2ip . then the Sawada's eigenmode method [131 E] is useful. 
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6 Basic equations for determination of parameters k 
and / in Lagrange multiplier terms 

The constraint term H' (13. 9p should vanish in the physical fermion space but its classical 
part in Eq in (15.21) might not. Thus we have an unsolved problem to determine the unknown 
parameters k a p and l a p in M a . They cannot be determined in the classical limit only. The 
determination of them requires taking the quantum mechanical fluctuations into account 
but has a difficult prescription. Instead, here we attempt another approach which has been 
proposed by one of the present author (S.N.) Ref. [9]. Following Eq. (6.1) in Ref. [9], using 
the relations 

c a % (G)=0, ct$oo(G)=ri$oo(G), ®* 00 (G) = <0 \U(G)\ 0>, 4 = -<, (6.1) 
we obtain the identities 

Ca U(Q) | > = {-r a +r a r^c\- q< c\) ■ U(G) | >, c\ U{Q) | > = - C t • U(G) | >, (6.2) 

which are just the identities given by Eq. (6.2) in Ref. [9] and whose proofs are given in 
App. A where we have used U{Q) \ 0>—U(G) \ 0> and (12.71) . repeatedly. It is shown that on 
the U(Q) \ 0> the fermion operators c a and satisfy exactly the anticommutation relations 

(clcp+cpct)U(g)\ 0> = 6 a p-U(G)\0>, ] 

/ \ ) (6-3) 

(cl l cl+cy a )u(g)\o> = (c a cp+c l3 c a )u(g)\o>=o, J 

for details of the corresponding proofs see App. A. The relations (I6.3p . however, play no 
role to determine the parameters k a p and l a p. At a region very near z=l where we have 

/ 1-1-2 

no effects of the unpaired modes discussed in the previous sections, y — — (= f) is approxi- 
mated to - — ^-j (0<g£~l). Then introducing the approximate function for the vacuum 
function which has first been proposed in Ref. [9], 

n (G)^{^U9)+mo(G)}^+^n o (G), $* M (g) = <0\U(g)\0> = ^* m (G), (6.4) 
we adopt an approximate SO(2N+l) WF expressed as 

U(G) \0> = % (G) (1 +r a ci ) exp {^q^M | > = U(G) | >, 

U(G) | > = U(G) | Q > = ~ U+ - ) U(G) 1 > . (equivalent to Eq. (6.6) in Ref. [9]) 

^oo (G) 2 \ pj 

We summarize some results in Ref. [9]. Using (16. 5 p and ( lA.6p . we get approximate identities 
c a U(Q) | > « c a C/(G) I > = ~ L*+ ~) c a C/(G) | > 



^(^-^)F(i+^)" 1 +^(i+^)- 1 } Q ^(G')|o>, 



4C/(S) I > « ct U(G) I o> 4^(G) I o > 



►(6.6) 



+ l(^)W 1+(?(?t ) ~ 1 -r T q j a+qq j r 1 } a u(G)\o>. 



With the aid of (1A.8I) . successive operations of the differential form (lA.lj) for the fermion 
operators c Q and on the approximate identities, though the formulas for such successive 
operations in Ref. [9] have been incorrect, leads to the correct formulas 



clcpU(G)\0> 



cp4,U(G)\0>: 



44c/(G)|0> 



1 



1 



-\c\ Cfi +-[i--\[l 



1 



1 



2 I <W+ 



1 



2^- 



32fV £ 

Z7(G?)|0>, 
3 



2 J "« 



32f 



Mt~, ( 









4) 



a /3 oo 7,4 



32f 



+ 



-7 ) (<cj + xjct 



- "7 I X* xl 



U(G)\0>, 



(6.7) 



C/(G)|0>, 



c a c p U(G)\0> 



1 

32/ 
1 



3 > r 

. I ry ry 



C/(G)|0>. 



Using the set of (16 .7p and (I6.8p . operations of the anticommutators of the fermion operators 
c Q and on the approximate U(G) \ 0> of the U{Q) \ 0> are made as follows: 



( C t C/3 + C/3C t)t/(G)|0>: 

7. 



i 



+ T 

■ 1 x a x /3 



2/ 



32/ 4 

c a cp + c p c a )U{G) \ 0> 



Fp~^ I ry i /'•"I" 

J, o t -/3 TX /3 L a 



^(C)|0>, 



E7(G)|0>, 



Z7(G)|0>. 



(6.9) 



We should demand that the expectation values of the anticommutators between c a and c^ a 
by the approximate U(G) \ 0> of U(Q) \ 0> satifsy the following relations: 

<0\W(G)(cicp + cpci)U(G)\0> = 8 a/Sj 
<0\W(G) ( C t c t +c t c t^( G )| O > = <O\W(G)(c a c +c c a )U(G)\O> = O : 
from which and (lA.llj) . we can obtain the following approximate relations: 



(6.10) 



<0\tf(G)(cicp + Cpci)U(G)\0> = ~y + ^ 



1 



1- 



2^ 



2 I <W " 



4^ 



(xj< C a >* G + X* a < Cp > G ) 



16f 



(6.11) 
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<0|C/t(G)(44+4ct) U(G)\0> = - 



32f 



<0\W(G)(c a cp+cpc a )U(G)\0> = ~ 
~^f(f~^) ( x * <c P >G+x J <Ca>G ) 



«0. 

'I 

>0, 



32f 



> (6.12) 



which are rewritten in compact forms as 

(3^ 4 - 2f + 1) 5 a/3 + — 2 ^ 2 - 3) x*^ = - ^ 2 - 1) {xp< c a >* G +<< ^ > G ) 



2^ 



^ — — J(a;*< C/3 >g+^^< c Q >g) and the complex conjugate. 



>(6.13) 



We call the demand mentioned above a quasi-anticommutation-relation approximation for c Q 
and c£, proposed first in Ref. The quasi- anticommutation relation- approximation may 
work well at the region very near 2 = 1, i.e., f = 1. To see more clearly the relation between 
SCF parameter M and Lagrange parameters k and Z, we give another form of SCF condition 
for M related directly to k and I. For instance, multiplying by kp 1 and lp^ the first equation 
of (I6.13P and its complex conjugation, respectively and summing up with respect to 7 and 
using the definition of M a in ( 13.12p . the conditions to determine the unknown parameters 
in the Lagrange multipliers for f 1 are given as follows: 



3/ 4 -2^ 2 -l „ . _ 1 f-3 



f-l 



(k 0a x a + l l3a xl) J r 



4f f-l 



•K<x I ( ktfj — I p-:X^ ) 

2 : 



x*<c a > G ) X] 7 ^7 x 7 _ ( x a< c a > G) J2- y kp- f x J = \x a \ 2 Mp (not summed for a 
1 



3" 



► (6.14) 



y^ 7 {kp^x-y+ip^x*^ 



- (x* a <c a > G ) Y.- y k (3'r x 7~( x »<c a >G) Y j1 1 M x * = \ x o\ 2m P (not summed for a) 

(x^<c^> G ) E7^7 X 7-( x /3 <c /3 > G)E 7 ^7 a; 7 = l a; /3| 2M a ( not summed for /3), 

3 N 



- 2^ 



7 

-1 



> (6.15) 



- {x* l3 <cp> G )Y J1 k ai x 1 -(xp<cp> k G )Y J1 l a 1 x* = \xp\ 2 M a (not summed for /3), ^ 

and their complex conjugations. To write the conditions in symmetric forms with respect to 
a and j3, we add (16. 14ft and ( 16.151) . Then for f^l we have 

2 ^ 2\ - ■{feA+'g/J + fe+^)} 



1 ^ 2 -3 



4^ 2 / 2 -l 



"cry "^7 



"0:7 A 7 



E 7 (^7^7 + ^7) + I X P 1 2 ^7 ( ^7 X 7 + ^ 
-X a < C a >* G YlfiPl X -r- X a< c a>G E 7 W^* - X ' P< C P >*G Y^4*aj X 7 ~ X *p< C P >G 

= \x a \ 2 Mp + \xp\ 2 M*, (not summed for a and (3). 
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(6.16) 



1 i 2 -3 



2/ 



2 ;2 



-❖^^M-vo^W-^^a^-v^tEM (6 - 17) 

= |x a | 2 M / 3 + |xa| 2 M a , (not summed for a and 

These conditions were first obtained in Ref. [9] but their expressions have been found to be 
wrong. Now we have their correct forms. Thus we could reach an ultimate goal to determine 
self-consistently the Lagrange parameters k a g and A/3- This may play crucial roles for a 
unified description of Bose- Fermi type collective excitations at a region very near ^ 2 =1, since 
they are the necessary and sufficient conditions to determine the unknown parameters in the 
Langrange multipliers. To gain trust in this approximation, we should compare it with a 
result of other Bose-Fermi theory in a concrete model. Toward a possible determination of 
them for f 1 m 1, it is better to make rearrangements of ( 16.1 6p and ( I6.17P with respect to k 
and / in the following forms: 
3^ 4 -2^ 2 -l 



■ (k/3 a x a -\- k a pXp) 



/ 1 i 2 — 3 \ /i^ 2 _3 \ 

+ (^2 I Xa 1 2 ~ X ° < ° a >ib ) ^ k ^ X " 1 + \4f f~l \ X P\ 2 - X l <C P>G) £7^7 

+2 3^-2^-1 (6.18) 

+ (^2 I ^ 1 2 - < < C « £7 l(h X y + (j^2 I ^ 1 2 - ^ < C /3 >fj) £7 ^ 7 ^7 



\x a \ 2 Mp + \xp\ 2 M*, (not summed for a and 0), 

j2 o \ / i j2 



Q 2 / 2 -l ^"l 2- ^^ Ca >g ) £7^7 - Q 2 ^ \xp\ 2 ~xl<cp >Gj Y,$cn* 
f 1 f-3 A f\ f-3 A 



= \x a \ 2 Mp + \xp\ 2 M a (not summed for a and (3), 

and their complex conjugations. Owing to the symmetric forms with respect to a and j3, 
(I6.18P and (16. 19[) can be separated into two parts. One is given as follows: 

/2 -i / 1 v2 



3g-2g-l / 1 ^-3 2 _ A 

f-1 \Aff-l / 7 

+2 3/_2^_J. .^ x * + ^^^| XQ | 2 _ x * <Ca>G ^^^ x * (6.20) 



= |x a | 2 M3, (not summed for a and /?), 

/ i ^ 2 _3 \ / J ^ 2 _3 \ 

~" V.2^ 2 ^ 2 1^°^ ~~ A< c a >G J ^7^/97^7 — y^ 2 "^ _ C a >g| E^^S (6 21) 

= |x a | 2 Mg (not summed for a and (3). 

The other is obtained by taking the complex conjugation and making the exchange between 
the indices a and /3 in ( 16.201) and only by making the exchange betwenn the indices a and 
(3 in (16.211) . respectively. 
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We are also able to consider another approximate SO(2N + 1) WF U(G)\0> at a 
region very near z = 0. At that region the SO(2N+l) WF has a peculiar feature 
that brings out the largest contributions from the unpaired modes owing to the relation 
x^Xq + z 2 = 1. We have never experienced such a physical situation. Then, an investiga- 
tion of such a interesting problem must be meaninigful. At the region very near z=0, the 



f, is approximated to 



WF expressed as 



V2 V 



flUz 



z~0). Introducing an approximate vacuum function 

r 



2f 



$oo we adopt an approximate SO(2N + l) 



Using (16.221) and ( 1A.6I) . we get approximate identities 



+ ^(C)|0 



>. 



c a U{G)\0>t 



V2\ U 



c a U(G)\0>- 



1 ( 1 



27! ^-^{^(^^^(^^^^(^I >> 



ciU(G)\ U>*—\f + -)cl U(G)\ 0>+:77 =^--){rt(l+ggt r i_ r T g t (1 ^ g t r i} Qf/(G) | >. 



2^ 2^ 



(6.22) 



>(6.23) 



With the aid of (1A.8I) . successive operations of the differential form (1 A. 1 1) for the fermion 
operators c a and 4 on the approximate identities lead to the formulas 



clcpU(G)\0> 



c p c{U(G)\0> 



clclU(G)\0> 



1 

7! 



i 



v/2 
1 



1 



2fi 



I6y/2f 



3 

2/. 



4^2> 2 



1 

2^. 



1 



1 



T t * 
1 







2^ 



1 



f/(G)|0>, 
3 



v^( /+ 2/. 



4V2^ 2 



1 

2£ 



^16^ 

C/(G)|0> 



★ T 



* T 

_ %KJ f) 



H6.24) 



+ 



16 v^/ 
1 



3 

2£ 



4V2~f 



2jt 



c a c 9 U(G)\0> 



16\/2> 4 
1 



2/. 



T 

XqiX p 



4V2~f 



Z7(GQ|0>, 



C/(G)|0>. 



(6.25) 



Using the set of (I6.24p and (I6.25p . operations of the anticommutators of the fermion opera- 
tors c a and 4 on the approximated U(G)\ 0> are made as follows: 



(tiep+epti)U{G)\ 0>' 



8VW 2 ^ 



* T 

/ >"» n rr> 



44+44Wg)|o> 



2^' 



4^ 



16^2/ V r 2^ 
and the complex conjugate. 



f 2^ ^ ^ 7 2^ 



2^ 



U(G)\0>, 



U(G)\0> 



(6.26) 
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The expectation values of the anticommutators between c a and c a by U{G) \ 0> should sat- 
ifsy ( 16.1 Op . from which and ( 1A.11I) . we can obtain the following approximate relations: 



<0\W(G){cl c p + c p cl)U{G)\ 0> = ^=U+ 



V2X 



W2U 



2/ 



1- 



2^ 



a/3" 



(xj< C a >* G +X* a < Cp > G ) 



8y/2f 



3 

2£ 



<0|m( G )(ct c t +c t c g^ G) |o> = _^ + -j^- w ^- 
1 



2 £ 



>(6.27) 



< C/3 >* G +^fl< C a >G 



:0 and the complex conjugate, 



which are rewritten in compact forms as 

2 (6^ 4 - 3f + 1) 6 aP + —2 (2^ 2 - 3) x* x„ = (2^ 2 - 1) (x p < c a >* G +a£ < c p > G ) , 

1 

2 (2^ 2 — 3)x*x^ = (2^ 2 — l) (x*<cg>G +x^<c q >g) and the complex conjugate, 

8^ 

which is converted into the following conditions to determine the parameters k and I: 



(6.28) 



,6^-3^ + 1 



1 2/ 2 -3 
Af 2f-l 



x n 



E 7 (^/37 X 7 + ^7 X 7) 



2^-1 

- (Xu<c a > G )J2^hi x *-( x a <c a>h)Y,^ k l3i x i = \ x a\ 2 Mf] (not summed for a) 



1 2^ 2 



> (6.29) 



x r 



8f2f-l 

- {x* a <c a > G )Y, 1 kp 1 x^-(x a <c a >* G )Y, 1 hi x * 1 = \ x a\ 2M p (not summed for a) 



,6^-3^ + 1 
2^ 2 -l 



1 2f 



Af 2f-l 



\ X p\ 2 E 7 {ka~/ X *+la"/ X -y) 



E 7 fc a 7 < - ( X P< C/3 >g) E 7 * 



7 " a-y X "/ 



\xp\ 2 M* (not summed for /3), 



( 

1 2^ 2 — 3 

Si 2 2^ 2 — 1 E 7 (^a7 X 7~'~^7- r 7 ) 

- (a^<C/3> G ) E 7 fc a7 a; 7 -(x / g<c^>G)Z] 7 ^7a; 7 =|a; ( g| 2 M a (not summed for /?), 



> (6.30) 



and their complex conjugations. Parallel to the way to get the conditions ( I6.16P and (I6.17p . 
from (I6.29p and ( I6.30p . we also have the symmetric conditions with respect to a and j3. 
Toward a possible determination of the unknown parameters for ^ 2 ^ |, we make again re- 
arrangements of the conditions with respect to k and I in the following forms: 



+ 



2f 
1 2^ 2 -3 
4f 2f-l 



6i 4 -3i 2 + l 



+ 



2^-1 
1 2^ 2 -3 



(6.31) 



a<C Q >G E 7 W< + 



1 2^ 2 



Af 2f-l 



|x^| 2 -x /3 <c / 3>* G )E 7 ^ 7 



2 /1/T* 
a ; 



(not summed for a and /J), 
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(6.32) 



= \x a \ 2 Mp + \xp\ 2 M a (not summed for a and j3), 

which can also be separated into two parts. One is given as follows: 

y4 n j2 , 1 / 1 ov.2 



6f-3£ + l ( 1 2^-3 , 2 _ A 

6^ 4 -3f + l . / 1 2^ 2 -3. I2 . \_ , . (6.33) 

+ 2 2 ^2_ 1 W + (^2 [ K I -< < Ca >gJE 7 ^ 

= |x Q | 2 Mg, (not summed for a and j3), 



_ f 1 2 r~ 3 1 i2_ * Vi- / 1 2f -3 \ ^ 

W 2 2/-l X « <C « >G )^^ \8f2f-l lXal x - <Ca>G )W^ x i (6 34) 

= |x a | 2 Mg (not summed for a and (3). 

The other is also obtained by the same way as the previous one. We are now in a stage to 
determine the unknown parameters k a @ and l a p in the Lagrange multipliers under the quasi 
anti-commutation relation approximation for the fermions. In the following section, as an 
illustrative example, we will treat the superconducting toy-model used in the previous Sec. 
4. The determination will be a great step to the completion of the present method. 
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7 Solutions for unknown parameters k and / 

First we give solutions at the region very near / 2 = 1. Using the Hartree eigenstate w a i, 
EHB eigenvectors U{ and Vi and the solutions for the unpaired-mode amplitudes Xi and Xj, 
which are already obtained in the previous Sec. 4, (I6.20p and ( I6.2ip are rewritten as 
3i 4 -2i 2 -l , _ + 3i 4 -2i 2 -l , _ + 

2 — -^Zi — kfSa ^ XiW *i +2 — — Pa ^ XiWai 

+E iY {(^2 - \ + M 2 ) x *' x ? ~ \i u >v +u i >v*,)x*a)jw ai/ w* aj ,J2~ ( hi E* x i w *i 

;.. (not summed for a and 0), 



(If 


^-1 2 
f 2 -3 1 




f 2 -l 2 








f-3 1 


r 1 * 


/ 2 -l 2 

; 2 -3 i 



(7.1) 



H x i%i 



]( 1 ^ 3 1 , i \2\ * 1 r * , *w 1 * / * ( 7 - 2 ) 

= J2ij w aiWaj Ei M i w %- ( not summed for a and 

Subtracting ( 17.2ft from ( 17. ip . we can get a simple relation 

3/ 4 -2/ 2 -l 



x (E 7 VyEi x i W 7i + E 7 ^7Ei ^*w 7 jj = 0, (not summed for a and (3), 

or adding (17.11) to (17.21) . we can reach another result 
3i 4 -2i 2 -l 

2 -72— \ (^oEi ^<t + ^cEi S >m) 



1 i 2 -3 



E iy '4^' w «'^v(E 7 ^7Ei^^+E 7 ^7Ei4 w 7i) 



4^ ^ — 1 J ^ V 7 ' l ^7^'^* * '7 (7 4) 

+ Ei'j' jd^'l 2- K'P) x i' x r - \ (u^Vif+Ui'V*,) X*,Xj> + ^ (upVf+TMjtVp) X^Xj^Wai'W^, 

x (^y k ^Y J i x i w *i-Y, 1 hiY J i x i w i^= 2 Y,i J w ( not summed for a and (3). 

The other result is obtained by taking the complex conjugation and making the exchange 
between the indices a and 0. 

Next we give solutions at the region very near f=\. (16\33|) and (l6\34j) are rewritten as 



2 ' 

% X i W ai 



6^-3^+1 6/ 4 -3^ + l 

— 2/ 2 -l ' /3q ^ ai — 2/ 2 -l — pa ^ 

f/ i 2i 2 -3 1 \ 1 1 

+E»'i|^ 4 ^2 2^_^ ~ \ + 1 v i' I 2 ) ~ ^ w ai '< 3 -E 7 ^7Ei s >7i 

= J2ij w aiWaj Ei ^i^ i3 (not summed for a and /3), 
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EiyjY 

2i 2 - 1 ~ \ + ^) X * i,Xj ' ~ 5 + ^*0^ ,C J , } W «* ,W ^7 Z ^< x * w r ^'^ 



1 2i 2 — 3 1 \ 1 

2 



\f2f-l 2 

= J2ij w aiWaj Ei Mi^i) ( not summed for a and 

Subtracting (17.6jl from ( 17.51) . we can also get a simple relation 



2^-1 



f/3 2i — 3 \ 1 1 1 

x (E 7 k M Ei X ^ 7 i + E 7 hi Ei s i w 7«) = °j ( not summed for a and (3) 

or addting (17. 5p to (17. 6p . we can also reach another result 
6/ 4 -3^ 2 + l 



og' 



2f-l 



■ ( fc /3«Ei ^<i + WEi X i W <*i) 



+ 8^2^T Eiy^'^'w*'^ (E 7 ^7Ei^^+E 7 ^7Ei^>7 V ^ 7 g ^ 

+ Ei'j'jd^'l 2 -!^'! 2 ) x *' x f -^(ut' V i'+ U i' V i') X i' X j' + ^( U j' V j' + U j' V j') X V X 3^ W oi>'Kg> 

x (E 7 fc /37Ei^Ki _ E 7 ^7Ei^ w 7*) =2 Eij w Qi w aiEi^^i ) ( not summed for a and 0). 

The other is also obtained by the same way as the above one. 

First we treat the region very near ^ 2 = 1. Subtracting and adding complex conjugation 
of (J73D and (OD from ((£3) and (J72J and to fl£3) and 1[73)» . respectively, we have 

3/ 4 -2^ 2 -l 



~^Ei'j 

X 



2 { (^«T^a)Ei^<+ ftaa^jEi x >cci) 

E 7 (^7T^ 7 )Ei^ W 7* + E 7 (^7 = F^ 7 )Ei :r i u; 7i} = ° ; ( not summed for a and 
Eiy j(- ^2 ^2— J + 5 - \ v f 1 2 ) 4^'+^ (upVf+UjiV^x^x^jw^w^, E 7 ( fc /3^ 7 )Ei ^7* 

+Eiy|^- A ^2 J| +^ - M 2 ) 4'^'+^ w m' w ^E 7 (^^ 7 )Ei^>7» ^ 7 ' 10 ^ 

= ^ijWaaWajY^iiMiWp^MiWpi) , (not summed for a and 

To solve (j7.9p and (I7.1(jp . we make the following bold assumptions: 

E 7 Ei x iK = E, M iW % E 7 ihi^ k h) E x iK = ^ E, M* WfH , 

(7.11) 

(not summed for /3), 
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E< 



1 ^ 2 -3 1 



„ 17 1 f-3 1 . l2 \ 1 . + ^ * 1 * (7.12) 



= - Y,ij w aiU>aji (not summed for a), 

in which the second assumption (I7.12p is proved to hold if the Hartree state is a real eigen 
state. We, however, use this assumption even if the Hartree state is a complex eigen state. 
Substituting (17. lip and (I7.12p into ( 17. 9p . we reach the following equation: 

y4 r, ;2 



1 " { {k/3aT l} a ) Ei %i<i + {ha T k*p a ) Ei x\w ai } 



(7.13) 



f 1 i -3 1 

= Ej U 2 + 4?^ ;E ^J WaiW * a ^^ MiW ^ M ^ W ^ ' ^ not summed for a and 
Equating each term in the l.h.s. and the r.h.s. of (I7.13p . respectively, then we have 

(not summed for a and /3), 
in which making summation with respect to a and using (17.111) . (I7.14p can be converted to 



■ 3 ^ 2 -^Ei2 + * fe^k, fcW^I. (7.15) 



from which we can obtain the relation between the original variable z and the number of 
single-particle state N as 

z 'L +1 _ (2W)i p^L 1 _i i (7 . 16) 

fortunately, which is a little bit smaller than 1 for large N. Further we study another case. 

Contrasting to the above case, subtraction and addition of complex conjugation of (17. 8p 
from (17. 8p and to (17. 8p . lead to a very simple equation 

3^ 4 -2/ 2 -l 



2 TTZi ^ ^ T ^«) ^ Xi ^ + ^ Tk lo)T,i %>ai } 



= 2 Y^ij w oLiWajY^i{ M i w %^ M i w Pi) 7 ( not summed for a and (3). 

In which we have used the trivial relation derived from ( 17. 12ft . Substituting (17.111) and (17.121) 
into ( I7.17p . we obtain 

y4 n i2 



2 ^ ft I 1 ' { ( k ^ l h) Ei X i<i + (ha T ^a) Ei X >ai } 



= 2 ^ijWaiWajYtiiMiWpiTMfwpi) , (not summed for a and /3), 

which is just identical with (I7.13p . 
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Next we treat the region very near f 2 = \- Subtracting and adding complex conjugation 
of (17. 7p and (17. 6p from (17. 1\ and (17. 6p and to (17. 7p and (17. 6p . respectively, we have 

x {E 7 ( fc /37=F^ 7 )Ei :E i W 7i + ^7(^7 = F^ 7 )Ei a; >7*} = °> ( not summed for a and /?), 
[( 1 2i 2 — 3 1 \ 1 1 

' 1 2i 2 — 3 1 \ 1 1 (7 20) 



= Y.ij w aiWaj^i( M i w %TM*w )3i ) , (not summed for a and /?). 

Instead of the second assumption (17.121) . we here make the third assumption 
{( 1 2i 2 — 3 1 \ 1 

„ 17 1 2^ 2 -3 1 . .\ . 1 . . . 1 . (7.21) 

= ^i'i ^£i 2 f-l ~ 2 + \ Vi y Xi ' Xj '~2 ( u i' Vi ' +u i' v i') x i' x ^ w o-vw ajl 

= ~ J2ij w aiWaj, (not summed for a). 

Substituting (17.111) into (17.191) . we obtain 

2 2i 2 -l {i k ^ l h)^i x *<i k} a ) 2~2i x*w ai } 



(7.22) 



r i 2i — 3 1 

= Y ~ ^2 2 ^2 1 X t x jj w aiW* aj J2i (MiWpiTMfwpi) , (not summed for a and (3). 
Equating each term in the l.h.s. and the r.h.s. of (I7.22p . respectively, then we have 

2 & \f-\ l ' fe^JEi XiW *i = fe-^7y^ x i x ^ w * w Z^ M i w h> (7.23) 

(not summed for a and (3), 
in which making summation with respect to a and using (17.111) . (I7.23P can be converted to 

2 6£ = 3£+l I 2 _J_^U 5 .. (7 24) 

2^ 2 — 1 8^ 2 -l 4 J J tv [ ' 

from which we acquire the relation similar to (17. 16[) as 



1 / 264 fi 

*=— 48iV-15-(8iV-15) 4 /l , (7.25) 

22 \ 1 ; y (8iV-15) 2 J 8JV-15' 1 } 

whose value is very near zero for large N. We study another case. Subtraction and addition 
of complex conjugation of (17. 8p from (17. 8p and to (17. 8p . lead to a very simple equation 

2 ^ 2f-\ l ' ^ kpaT ^ ^ XiW * ai + ( l(3a T k P°) ^ X * Wai } 

+ 87 2 2^1 ^iT x i' x i' w ™'K'&( k ^ l h) Ei^*i + E 7 (^ 7 : F^ T ) E<^>7*} ^ ? ' 26 ^ 
= 2J2ij w aiWajJ2i{ M i w hT M i w l3i) > ( not summed for a and (5). 
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In which we have used the trivial relation which is also derived from (17.211) . Substituting 
(17. lip into ( I7.17p . we obtain the equation equivalent to (I7.22p . This situation is quitely the 
same as the one of the previous case. Using (17.131) and (I7.22p . the Lagrange multipliers k a p 
and l a /3 are determined in the following forms in which the indices a and are not summed: 

fea^ a =ciEij( 2 +C2^^i) Wf,iW* p S MlW f , k fio ^l=c^ i3 i2M: 2 x\x j ) w m < ^ MlW } 1 , (7. 27) 

Z^i X i W pi Z^i X i W ai 

first of which is brought by exchange of the indices a and (3 of the second which corresponds 
to (I7.13P and (I7.22p . The complex conjugation and exchange of the indices of (I7.27P bring 

k a ^p a =c 1 J2 ij {2+c 2 x t x j)w ai w* aj -^ — r~^' kf 3o ?flf 3a =c^ i3 (2+C2X*x j )wpiWl j ^4 J„ Wm .(7.28) 



Y,i x >pi 



We have used the properties kR a =k^ g and l/3 a =lai3- The coefficients C\ and c 2 are given as 



ci = 8- 



21 1 
~2 2N+V 



( 



c 2 ■ 



7 



(2N+1) 



1- 



7 



1 



V 



8 2iV+l-^ 



8 8iV 2 -21iV+18 



Cl = 3" 



c 2 ; 



1 (8JV-15)(4iV-9) 



(7.29) 



8iV-15 ' ^ 6 8iV-9 ' vr 2y 
We here should notice that by only the equations ( 17.271) and ( 17.281) the unknown parameters 
k a/ 3 and l a p cannot been determined. Such a determination is possible if the equations (I7.27P 
and (17.281) are subtracted or added with each other. Really, subtraction the second from the 
first in (17271) and (1728) leads to 



k a/3 -kp a = C 1 J2M+C2 x t x j)w^iWi ^ M%W ™ -C^2J2+C 2 x i x j)WaiW* ^' 



Ei x i< 



k a f} kp a : 



J Ei X i W ai 



respectively. From these relations, we can determine the unknown parameter k a g as 



E.M^ l 

( 

2 



i X i W c 



On the other hand, adding the second to the first in (17.281) yields 



k a p+kp p^l Ba =c^ b (2+c 2 x*x j )w ai wlp=± — ^-^+cEi/ 2 + c 2^*^i)w/3i^- 



from which we are able to determine the unknown parameter l a g as 



Ei x >fn 



(7.30) 
(7.31) 

(7.32) 

(7.33) 



Tlp a =^^2ii 2 +c 2 x*x j )w ai w 



E, Mtvpt , 3 



Ei Z>m 4 



t-l c iEi ? (2+c 2 x*a; i )^W ( 



E, Mtw a 

Ei 



1 V- (n r * \ * Ei 



1 E M w* ■ 

r c iE t /2+c 2 x*x j )w / 3iW^. ^ - -fg = T^- 



(7.34) 



To verify the validity of k a p ( 17. 32ft . we make the explicit form of A^q, from fc aj g by exchanging 
the indices. We calculate k a p — kp a and to the result of which further substitute (17.271) and 
(I7.28p . Then we can get k a/ 3 — k 8a again. This means that the expression for k a g is exactly 
valid and does not cause any inconsistency. By the direct input of the explict expressions 
(I7.32p and (17.341) into the SCF parameter M a difined by (13. lip , we could reach our ultimate 
goal to determine self-cosistently the unknown Lagrange parameters k a p and l a p. 

Finally summarizing shortly, in the present manner we could reach considerably good 
solutions for all unknown parameters k and I at both the regions very near ^ 2 =1 and ^ 2 



2" 



We will apply the present approximation method to concrte problems which should be solved. 
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8 Concluding remarks and further perspectives 



In the present new description, the extended HB (EHB) theory for a fermion system with 
N single-particle states has been derived from the extended TDHB (ETDHB) theory. The 
EHB and ETDHB theories have been constructed from a group theoretical deduction start- 
ing with the fact that the fermion annihilation-creation and pair operators form a Lie algebra 
of the SO(2N+l) group. The induced representations of the SO(2N+l) group was found by 
means of a group extension of the SO(2N) Bogoliubov transformation for the fermions to 
the SO(2N+l) tranformation group. Embedding the SO(2N+l) group into the SO(2N+2) 
group and using the boson images of the SO(2N+2) Lie operators, we have obtained the 
ETDHB equation from the Heisenberg equation of motion for the boson operators. We have 
expressed its final form through the variables as the representatives of the paired modes and 
the unpaired modes. From the ETDHB theory we have derived the static EHB theory in 
which the paired modes and the unpaired modes are treated in an equal manner. The EHB 
theory applicable to both even and odd fermion systems is also the SCF theory with the 
same capacity to provide a mean field approximation as the usual HB theory for the even 
fermion systems. We have obtained a new solution with unpaired-mode effects. However, 
it includes the unknown parameters which originate from the Lagrange multipliers involved 
in the image of the Hamiltonian in order to select out the physical spinor subspace. They 
cannot be determined in the classical limit only, and a complete determination of them re- 
quires that the quantum mechanical fluctuations are taken into account. We have, instead, 
determined the parameters with the aid of the quasi- anticommutation relation approxima- 
tion for the differential form of the fermion operators. Both the EHB eigenvalue equation 
and quasi- anticommutation relation approximation provide a group theoretically transpar- 
ent kinematical frame, which is able to describe both the paired and unpaired modes in 
many fermion dynamics. They work well in both even and odd fermion systems with strong 
collective correlations in which the effect of the unpaired modes is remarkable. They play 
important and crucial roles for a unified SCF description of the Bose-Fermi type collective 
excitations in several fermion systems, e.g., quantum dots [33]. We stress again that it is a 
very exciting problem to consider an approximate SO(2N+l) wave function very near z=0. 
We have never experienced such a physical situation in which the largest contributions from 
the unpaired-mode amplitudes occur owing to the constraint XaX a +z 2 =l and consequently 
most physical contributions arise from the unpaired-mode effects. This situation is in strik- 
ing contrast to the situation in which we have no effects due to the unpaired modes very near 
z—1. A study of such a problem can be expected to open quite a new field for exploration of 
an exotic fermion dynamics accompanying the stronger unpaired-mode effects. Particularly, 
in a forthcoming paper, we will investigate a concrete problem for such an object using the 
present superconducting toy-model. Then we will give detailed processes to the EHB eigen- 
value equation under the quas- anticommutation relation approximation very near z=0 and 
make analysess of solutions obtained by calculations. 

Hitherto we have no effective method to describe such Bose-Fermi type collective motions 
except the previous theory proposed by Fukutome and one of the present authors (S.N.) Ref. 
[31] in which the SO(2N+l) density matrix plays a crucial role. To study the interrelation 
between the previous theory and the present theory will become a very interesting project. 

Finally we will develop a group theoretical approach to the formation of the Lax pair of 
the SO(2N+2) top according to the ideas of Olshanetsky and Perelomov and Reyman and 
Shansky [20]- [22] in the near future. 
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Appendix 

A Differential SO(2N+2) Lie operators on coset space 
and operation of differential annihilation-creation 
operators on SO(2N+l) wave function U(Q)\0> 

Owing to <0|f/(G)|0> = <0|f/(£)|0>, we have | G> = U{Q)\ 0> = \Q >{Q <ESO(2N+2)). 
Following Refs. [23] and [21], expressions for c a and 4 are given in terms of the variables 
q a p and r a : q q q q q 

+r,^^ + ir a r i -q < )— — — +ir a — , = (A.l) 



Using (lA.ip we prove the identities ( 16. 2ft as follows: 
Ca U(G) | > = c a [/ (G) I > = c a |$* (G) (1 +r p cl)eM4 | > 







-r a U{G) \ 0>+^(G)l(r a r^-q o ^)c^e^ s< ^ c s\ 0>-(l+r /3 c^)q a (.r ri -^—e^ q -> sct -> c ^\0> 



i (A.2) 



r a C/ (G) | Q> + \ [ r a r^c\{l+r v c\-r n c\ 1 )-q < c ] ll +r v c\ -r,cj) | $£ (G)e3 W=s | > 

-g a ^4ct$5 (G)e5^ c M|0> 

r a t/(G)| 0>+r a r 5 4c/(G)| 0>-g aC 4C7(G)| 0> 

-ra+rorfcj-gafcjj -C/(G)| 0>, 
where we have used rgclr^cjj = 0. Thus we have given the proof of the first identity of (16. 2D . 
Further using r^c^r^c^ = 0, the second identity is also derived as follows: 
ct U{G)\ > = - c* |$* (G)(1 +r^4)el^44 1 >] 







% Q {G)l^-{l+r^c\)e^ c ^\ 0> + (l+77?cJ)r£-^- e ^ c W|0> 



9 



<9g, 



(A.3) 



= - C t$* (G)(l + r € 4-r^r 5 44) e l<M44| > = -cj,- C/(G)|0>. 
By successive use of (16. 2p and (16. ip . we get the formulas below 
cicpU{G)\ 0>=(5 aP -r p cl t +^44-5^44) -f/(G)| 0>, 

Cp ci u{G)\ 0=^4-^44+5^4 4) .f/(G)| o>, 4c^)| o>=-44-c/(G)| o>, 



c a cpU{G)\ >=|g a/3 -g a5 g^44 _ (^ r / 3_ra ^)4 + ^ Q « r / 3_ra ^) r ^4 c i ) ' f/ ( G )l > 
The first and last equations of (1A.4|) are proved as follows: 

cic p U(g)\ 0> = {4 (-r^+r^4-g^4)+(-r^+r^4-^4)(-4)}-[/(G)|0> 



(A.4) 



<W -^4+^44 -<?/3?44 ) ■ ^(c) | o >, 



CaC^(g)|0>: 



(r a r e -^) 



d ' 



-r P +r^cl-q K c\ 



WA.5) 



+ {^rp+r p r^c\-qp^{-r a + r a r^c\-q ai :cfjyU{G)\ 0> 
q a p-q^qpr l c\cl-{q a ^rp-r a q K ) c\ + (q a£> r p -r a q K ) r r] c\c ] \ -U{G)\0>. 
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Using another form of i.e., ^ = {l+r^x r ) 2 (x = (l+99^) *)> we gi ye useful formulas to 
evaluate expectation values of the anticommutators of the operators c a and in differential 
forms by an approximate U(G)\0> of the SO(2N+l) WF U(G)\ 0>, i.e., fl63|) : 



1 x T 
c ^ = -2^{ rT ( 1 + ^) _1 + rt ^( 1 +^) _1 } a = { rT ® rt } a = -J i 



c t^= ^ {r^l + gg t ) _1 -rV(l + W t ) _1 } a = { rt © rT } ( 

The proof of the first equation in flA.6[) is given as follows: 
d 



x 



l+z 



(A.6) 



1 1 



5 



9r J 2 1+rt X rV l+rt X r«9r* ( 1 + r ^W>) 2 ^ (Xr)a ' 
(9 19 

9 1 9 1 

( rar^-gaf) ^- ^ = - -f (r a r ? -g a? ) — (l +r^X/3 7 ?" 7 ) = {r f X^ " r - (QX*r*)} a , 



> (A.7) 



9rc 



(9 



9 



rgaer ^ o(i+ gg t) fe { (1+ggt) " 1 }^ 9(1 ^f )fc ^ 



Gathering these results, the first equation is obtained. Using the relation = — c*, the last 
equation of ( lA.lj) . the second equation is also proved. With the use of the same simbols as the 
ones introduced in ( 1A.6I) . we also give the following useful formulas first presented in Ref. El: 



c Q {r T ffirH 



-2 1 



^{rter T } a =2^1-^.^ + 



* T 
2~- 



2fi 



■fiat 



* T 



(A.8) 



The first and second equations of (1A.8|) are proved as follows: 

A{ r T 0r t^ = j g ( 1 + g t g )-l}^ ; 



dr* 



d 



r £-qai)-K- {r 1 ©? 4 } = r Q (r T (l+g+g) x } -{g(l + g f g) x } 



<9 



{r T ©r t } 



a/3 



* T 







m+oVrn d(i+gtg) 7V 

+r * r t <7 . ? lU+ata)- 1 ! g(1+gtg) ^ 

+ € ^ S d(l+qH)^ {[l + qq) I* dqK 



^(l+ggt) r | g (l + g t g )-l| + | r T( 1 + g t g )-l +r t g (l + g t g )-| a | r t g (l + g t g )-l| 



> (A.9) 



d 



<9 



-1^ v r^-^- {{l + q ] q) 1 }^-?af^r|r|g e a- 



( {r T (l + gtg)-l}^ + r t(l + ggt) r | g (l + g t g )-l}^ 

+ {r T (l + g t g)~ 1 +r t g(l + g t g)- 1 } Q {r T (l + g t g)~ 1 } /3 . 
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On the other hand, we have 
d 



d 



- (r* r\ - q^ ) — { r T ©r t } p = -r* { r\l + gt } p - ^+ { (1 + gtg) - 1 }^ 



_r £7r— { rT © rt }« = - r £ r JTT7-; — — 1(1" 



-i} + 



77 



^ dq. 



- t (l+gg t )r{(l+ g t g )-i} { r t(i +gg t)-i_ r T ? t(i +?? t)-i} a { r .T( 1+g t ? )-i} 



MA.10) 



9 



d 



/3 H <*i V 



d(l+qiq) 



. {(1+?t?) _ 1} aa^fe 



77 



+ Q^r*r^rlq eS 



d 



'9(l + gtg) 77 , IV- '-^ J5/3 ^ 

= r t (l + gg t )r-5 a/3 -r t (l + gg t )r{g(l + g t g)" 1 } a/g + rt {r t g(l + g t g)" 1 } /3 

+ {r T (l + g t g)" 1 g t -r t (l + g t g)" 1 } Q {r t g(l + g t g)" 1 } /3 

Gathering results ( 1A.9I) and ( 1A.10I) . the first and second equations of ( 1A.8|) are obtained. By 
using the relation c£ = — c* and the second equations of ( 1A.8I) . the last equation of ( 1A.8I) is 
also proved. Using the identities (16.21) and the two important relations < >g r% = and 

9a£<C£>(; = — 2 ( x +? a '*)a+ <c a the following expectation values of the annihilation and 
creation operators, c a and , are derived: 



<0|m(^) Ca f/(G)|0> = --^+-J|<c Q > G +^(x+gx\^--(^ + 7 )<c Q > G , 

<o|^) c t^(G)|o>=-i^+i > )< Ca >* G , (o«^i), 



KA.11) 



which have been used to compute the expectation values of the anticommutators between 
c a and ct by U(G)\0> (TOTj) . 
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